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ABSTRACT 

Aims. Using an AKARI multi- wavelength mid-infrared (IR) survey, we identify luminous starburst galaxies at z S 0.5 based on the 
PAH luminosity, and investigate the nature of these PAH-selected starbursts. 

Methods. An extragalactic survey with AKARI towards the north ecliptic pole (NEP), the NEP-Deep survey, is unique in terms of a 
comprehensive wavelength coverage from 2 to 24 /jm using all 9 photometric bands of the InfraRed Camera (IRC). This survey allows 
us to photometrically identify galaxies whose mid-IR emission is clearly dominated by PAHs. We propose a single colour selection 
method to identify such galaxies, using two mid-IR flux ratios at ll-to-7yum and 15-to-9/mi (PAH-to-continuum flux ratio in the 
rest-frame), which are useful to identify starburst galaxies at z ~ 0.5 and 1, respectively. We perform a fitting of the spectral energy 
distributions (SEDs) from optical to mid-IR wavelengths, using an evolutionary starburst model with a proper treatment of radiative 
transfer (SBURT), in order to investigate their nature. 

Results. The SBURT model reproduces observed optical-to-mid-IR SEDs of more than a half of PAH-selected galaxies. Based on the 
8 jim luminosity, we find ultra luminous infrared galaxies (ULIRGs) among PAH-selected galaxies. Their PAH luminosity is higher 
than local ULIRGs with a similar luminosity, and the PAH-to-total IR luminosity ratio is consistent with that of less luminous starburst 
galaxies. They are a unique galaxy population at high redshifts and we call these PAH-selected ULIRGs "PAH-luminous" galaxies. 
Although they are not as massive as submillimetre galaxies at z ~ 2, they have the stellar mass of > 3 x 10 10 M G and therefore 
moderately massive. 

Key words, infrared: galaxies - galaxies: starburst - galaxies: active - galaxies : evolution 

1 . Introduction which requires knowing the appropriate mid-IR spectral energy 

distribution (SED). In fact, sensitive IRS observations with the 
In recent studies of galaxy formation and evolution, the most Spkzer Space Telescope reveal that the mid . IR spectra of sub . 

massive and luminous galaxies have played a leading role. Most ^metre galaxies resembles that of over 2 ord ers of magnitude 

luminous galaxies, such as submillimetre galaxies and ultralu- lesfj luminous starburst galaxies , Valiante et al. 2007; Pope et al. 

minous infrared galaxies (ULIRG^ at high redshifts are now 2008; FaT rah et alJlTO . which have prominent emission fea- 

beheved to be the progenitors of massive spheroids, and high- mres of polycyclic hydrocarbons (PAHs). This raises 

light the early formation of massive galaxies (e.g. | Smailetal. | questions about the nature of most i uminous ga i axi es, and how 

|2004t | Chapman et al. || 2005t |Tacconi et al. || 2008|) . This supports they differ systema tically from low-z to high-z. 
the scena rio of anti-hierarchical or 'down-sizing' galaxy forma- 
tion (e.g. |Heavens!etld]|20^|C^ Multi-wavelength mid-IR surveys with AKARI provide a 

Mid-infrared (IR) surveys are a powerful tool for study- unique measure of the mid-IR properties of infr ared galaxies 

ing the most luminous galaxies at high redshifts, by pro- at high redshifts with a statistically valid sample. iTakagi et all 

viding the infrared luminosity function and the star forma- d2007al) demonstrate that AKARI/IRC all-band photometry is 

tion rate in a large cosmic volume as a function of red- capable of identifying the approximate spectral shape of the PAH 



shift (e.g IPerez-Gonzalez et al.1 l2005t iLe Floc'h et al.l l2005t emission, specifically the steep rise of flux at the blue side of the 

Babbedg eet all 120061; ICaputi et alj 120071: iGoto etal.1 1201 Oh . PAH 6.2 //m feature. Such a steep flux rise can only be accounted 

However, this work can suffer uncertainty from the /^-correction, for by PAH emission. Using this fact, we propose a new pho- 

tometric selection method for galaxies whose mid-IR emission 

1 Galaxies with a total infrared (8 - 1000 yum) luminosity of is dominated by PAH emission. This is possible if photomet- 

10 12 ~' 3 L Q . ric observations are sufficiently comprehensive enough to trace 
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the PAH emission feature in the SED, as in the AKARI/IRC 
all-band photometric survey. We call galaxies selected in this 
method PAH-selected galaxies. 

In this paper, we utilize the multi-wavelength mid-IR cov- 
erage of the NEP-Deep survey with AKARI, in order to iden- 
tify galaxies with strong PAH emission. Sections 2 and 3 de- 
scribe the data and sample selection, respectively. Our SED- 
fitting method using the optical to mid-IR wavelength range, 
and its results are presented in section 4. We discuss the prop- 
erties of PAH-selected galaxies in section 5 and summarise the 
results in section 6. Throughout this paper, a flat cosmology with 
Hq = 70km s _1 Mpc~' and Q.^ = 0.7 is used. All magnitudes are 
given in the AB system, unless otherwise stated. 



2. Data 

Here we use the multi- wavelengt h data set of the AK ARI NEP- 
Deep survey, described in detail in Wada et al. (2008). The NEP- 
Deep survey covers a circular area of 0.38 deg 2 , centered at 
a = 17 h 56 m 01 s and 6 = 66°33'48.4" (J2000.0), with all nine 
IRC bands, i.e. 2.4, 3.2, 4.1, 7.0, 9.0, 11, 15, 18, and 24 fim (N2, 
N3, N4, S7, S9W, 511, L15, L18W, and L24 in conventional 
band names, respectively). The 5 <x sensitivities in these bands 
are 9.6\ 7.5 1 , 5. 411 49, 58,71 117 121 and 275/Jy from N2 to 
L24, respective! v lWada et al.ll2008l) . Deep optical (BVRi'z' and 
NBlll) images with the Subaru/Suprime-cam (S-cam), reach- 
ing a limiting magnitude of B = 28 and z' = 26 AB mag, 
are available for a part of the NEP-Deep field with an area of 
27' x 34', i.e. one field-of-view of S-cam. We also have ground- 
based near-IR (J K S ) images with KPNO-2. 1 m/FLAMINGOS 
for the S-cam field dlmai etalj2007l) . reaching a Vega magnitude 
ofK s = 19.5. The ground-based near-IR images have higher spa- 
tial resolution (FWHM of 1.08") than that of IRC, and are there- 
fore quite useful to resolve source confusion in IRC images. In 
this paper, we concentrate on sources in the S-cam field only. 

We performed an optical spectroscopic survey of 242 mid- 
IR sources with R S= 24 mag in the S-cam field using 
Keck/DEIMOS (Takagi et al. in prep). Here we use the spec- 
troscopic redshifts for the calibration of photometric redshifts. 
The optical emission line diagnostics of mid-IR sources will be 
given elsewhere. 



3. Sample selection 

3.1. All-band-detected sources 

In order to evaluate and to maximize the new diagnostic ca- 
pabilities becoming available from the NEP-Deep survey, we 
have constructed a source catalogue of objects, detected in all 
nine IRC bands (hereafter all-band-detected sources). We have 
combined the 57, S9 W, 5 1 1-band images for source detection 
and used SExtractor dBertin & Arnoutsj 1 996|) to make an ini- 
tial catalogue. We then follow Takagi et al. ( 2007al) for photom- 
etry with IRC images, in which the aperture photometry with 2 
pixel and 3 pixel radii are used for near-IR {N2,N3,N4) and 
mid-IR images (S1,S9W,SU,L15,L18W,L24), respectively. 
Appr opriate aper t ure cor rections are then applied. We also fol- 
low iTakagi et al.l d2007al) for band merging. Starting from the 
centroid position of 5 1 1 sources, we searched for the centroid 
in the other images. If the angular distance between centroids 
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Fig. 1. A colour-colour plot of all-band-detected sources with 
15-to-9 //m and 1 l-to-7 //m flux ratios. Dashed lines indicate the 
colour cuts for selecting PAH-selected galaxies. Points with er- 
ror bars indicate the PAH-selected galaxies. 




Fig. 2. ll-to-7//m flux ratio as a function of redshift. Solid 
circles with error bars indicates the spectroscopic sample with 
5 7, 5 g,S ii > 30/idy and S 15 > IOOaJv. a) Flux ratio estimated 
from the SED template of (Polletta eTal] (120071) . The names of 
galaxies in the template are given in the legend, b) Same as a), 
but with the SBURT model with different starburst age and com- 
pactness of the starburst region (0=1.4, 1.6, 1.8, 2.0) , which 
controls the optical depth. The MW dust model is adopted. [See 
the electronic edition of the Journal for a colour version of this 
figure] 



1 Corrected from the values in Wad a et aTl d2008h . in which the 
adopted conversion factors for N2, N3, and N4 bands are not correct. 
The correct sensitivity is better by a factor of 1.47. 
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Fig. 3. Same as Figure^), but for 15-to-9/im flux ratio. 



is less than 3", t he objects are considered to be the same (see 
IWada et al.ll2008l Since we have not adopted a specific signal- 
to-noise ratio for the detection at all of the IRC bands other than 
S 11, the resulting catalogue includes less significant detections 
(< 3 cr) in some cases. Nevertheless, we include these sources 
in our all-band-detected sample, since the less significant detec- 
tions are rare, and the resulting SED is still well constrained. 

There are 1 100 all-band-detected sources in the entire NEP- 
Deep field, and 630 in the S-cam field. We searched for optical 
identification using sky positions in the N2 band, where the as- 
trometry is based on the Two Micron All Sky Survey (2MASS) 
and therefore most reliable. The resulting optical identifications 
were all visually inspected. We found that about 10 % of sources 
have ambiguous optical identifications. Excluding these sources, 
we obtained 568 all-band-detected sources in the S-cam field 
with unambiguous optical identification, of which 113 have op- 
tical spectra with DEIMOS. 

3.2. PAH-selected galaxies 

At the mid-IR wavelength region, emission from galaxies is usu- 
ally dominated by small and hot dust grains experiencing the 
temperature fluctuation. PAHs are believed to be the smallest of 
those, consisting of ~ 100 carbon atoms grouped into ~10 A in 
dimensions, and having emission bands. PAH emission features 
have sometimes been used as star formation indicat ors, which 
are prominent in sta rburst galaxies ()Brandl et al.l2006l) . but weak 
or absent in AGNs dWeedman et al.ll2005l) . 

Using the very steep flux rise on the blue side of the 6.2 /jm 
PAH feature, we attempt to identify galaxies whose mid-IR 
emission is dominated by PAHs, based on AKARI photomet- 
ric observations. The steep flux rise due to PAH emission can be 
characterized by the flux ratio at the rest-frame wavelengths of 
approximately 7 to 4 jj.m. For galaxies at z ~ 0.5 and 1, we could 
use 1 l-to-7 yum flux ratio and 15-to-9 /urn flux ratio, respectively, 
in order to identify galaxies with prominent PAH emission. For 
convenience, we hereafter use the term, the PAH-to-continuum 
flux ratio, for these flux ratios. Figure Q] shows a colour-colour 
plot of all-band-detected sources using these flux ratios. There 
are galaxies with the PAH-to-continuum flux ratio as high as 10 
for both ratios. We hereafter call galaxies with the flux ratio of 
> 8 as PAH-selected galaxies. 

For star-forming galaxies, it is likely that the flux at the rest- 
frame 4 //m is dominated by the stellar continuum, while PAH 
7.7 yum feature dominate the flux at the rest-frame 1 fim. The 



PAH luminosity is believed to be a good tracer of total infrared 
luminosity o f star- forming galaxies, and hence the SF R (e.g . 
iGenzeletal .11 19981; iRigopoulou et al.lll999b iFarrah et al.l 12007). 
Thus, these flux ratios are roughly proportional to the mass- 
normalized SFRs, i.e. specific SFR. We hence expect that PAH- 
selected galaxies are the best candidates of starburst galaxies. 

In Figure [2^, we show the 1 l-to-7 yum flux ratio as a function 
of redshift, using the SED t emplate of various galaxy types - 
the SWIRE template library (Poll etta et al.ll2007l) . This template 
confirms that the 1 l-to-7 yum flux ratio has a peak at z ~ 0.5, 
although no SEDs in the template can reproduce the flux ratios 
greater than 8. We note that the SWIRE template library has only 
a few SEDs of starburst-dominated galaxies. In order to explain 
the high flux ratio of PAH-selected galaxies, we may need more 
comprehensive SED templates of starburst galaxies. 

For such SED template of starburst galaxies, we adopt 
an evolutionary starbur st SED model with radiative trans- 
fer dTakagi et al.l l2003ai SBURT). This is an SED model for 
starburst regions, assumed to be spherical and have a con- 
centrated stellar distribution and homogeneous dust distribu- 
tion. SBURT adopts a ste llar population synthesis model of 
Kodama & Arimoto] ( 1 1997b as input. The amount of dust is cal- 
culated, based on the chemical evolution assuming gas infall. 
The time-scale of gas infall and of star formation is assumed 
to be O.lGyrB Thus, the SBURT model can be regarded as 
a population synthesis model with physically consistent treat- 
ment of dust extinction and re-emission. The optical depth of the 
system is controlled by the compactness (0) of the starburst re- 
gion. In the SBURT model, there are three choices of the extinc- 
tion curve, corresponding to the dust model for the Milky Way 
(MW), Large or Small Magellanic Clouds (LMC or SMC). The 
fraction of PAHs in dust is 5, 1 and 0.5 % in the MW, LMC, and 
SMC dust model, respectively. Thus, models with the MW dust 
has the largest PAH-to-continuum flux ratio. See iTakagi et al.l 
(l2003allbl) for more details of SBURT. 

In Figure |2j) and [3] we show the 1 l-to-7 /mi and 15-to-9yum 
flux ratios estimated from the SBURT model, respectively. We 
found that the MW dust model is required, in order to repro- 
duce high flux ratios of PAH-selected galaxies. According to 
the model, younger galaxies may have larger flux ratios. This 
is because the contribution from the stellar component to the 
4 yum continuum is less for younger starbursts. The flux ratios of 
greater than 8 are satisfied with the model younger than 0.4 Gyr. 
However, this is only one possible interpretation, since the dust 
model adopted has only three choices (MW, LMC or SMC), 
and limited capability of reproducing large variations of mid-IR 
properties. Specifically, this model does not take the ionization 
status of PAHs into accoun t, which significant ly affects the opti- 
cal constants of PAHs (e.g. lLi & D raine 200U). 

Among all-band-detected sources, we identify 56 (39) PAH- 
selected galaxies from 1 l-to-7 yum (15-to-9yum) flux ratio in the 
entire NEP-Deep field, of which 38 (18) lie in the S-cam field 
with unambiguous optical identification. We include sources 
with less significant detection at the rest-frame 4 fim, if the 2.5 cr 
upper limit is consistent with the selection criterion. In Table 
A. 1, we tabulate the IRC photometry of PAH-selected galaxies. 



2 The model SEP of starbursts scales w ith the time-scale of gas infall 
and star formation (Takaei et al. 2003a). Therefore, the starburst age 
from the SBURT model fitting should be considered only in the relative 
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Fig. 4. Colour-colour plot with N3-S 7 versus N2-N3. Solid cir- 
cles and small crosses indicate the good and bad-fit sample, re- 
spectively. Solid circles with error bars represent optically iden- 
tified AGN. 



4. Photometric redshifts and SED fits 

4.1. Method 

We fit the optical-to-MIR SEDs of all-band-detected sources 
with the SBURT model. We use the model with the starburst 
age of 0.01 - 0.6 Gyr, which is enough to represent the star- 
burst phase. The SBURT model ori ginally covers the wav elength 
range from UV to submillimetre. ITakagi et al.l (l2007bl) empir- 
ically added the radio component to reproduce the observed 
radio-IR correlation. This model does not include any emission 
from AGN. Thus, it is expected that the SBURT model under- 
predicts the mid-IR fluxes from AGN hosts, which us ually have 
hot d ust excess compared to pure starburst galaxies (Hu nt et al.l 
fl997h . 

We adopt the same SED fitting method as in ITakagi et all 
(l2007ah. usin g a standard x 2 minimization technique. 
Tak agi et all (I2007ah 



adopted relatively large minimum 
photometric errors, 20 - 30% depending on the band. Since 
the photome t ric c alibration of IRC is recently improved 
(Tanabeetal. 2008), we here set the minimum photometric 
errors to 15% for all IRC bands. As in ITakagi etal] d2007al) . 
we quadratically added an additional 20 % error for data at the 
rest-frame UV wavelengths, < 4000 A, in order to account for 
the uncertainty of the extinction curve. 

We apply the SBURT model for all-band-detected sources. 
By nature, all-band-detected sources are heterogeneous sample, 
even including stars. Since we use a starburst SED template, 
SBURT, for this sample, it is important to discriminate starburst 
candidates from other type of sources. To do this, we use the 
goodness of SBURT fitting and reject sources with no accept- 
able SBURT model in further analyses. We reject the best-fitting 
SBURT model if the value of resulting x 2 is so large that the 
corresponding probability is less than 1 %. 

4.2. Fitting results 

4.2.1. All-band-detected sources 

We performed SED-fitting for 568 all-band-detected sources 
using the SBURT model. We obtained good fits for ~ 40 % 
(211/568) of all-band-detected sources. On the other hand, us- 
ing the SWIRE template library, instead of SBURT, we obtained 
good fits for only ~ 10 % (69/568) of sources, although it covers a 




Fig. 5. Photometric versus spectroscopic redshifts. Solid circles 
indicate the redshifts of the good-fit sample, while small crosses 
are for the bad-fit sample. Dotted lines represent = ±0. 1 . 
Squares indicate the photometric redshifts of the good-fit sam- 
ple, based on the ground-based optical-NIR photometry only. 




Fig. 6. Correlation between N2 - AG colour and redshift for both 
photometric and spectroscopic samples. Open circles with error 
bars indicate the spectroscopic sample. Solid circles represent 
the photometric sample. Solid and dashed lines are the regres- 
sion line for photometric and spectroscopic sample, respectively. 



wide range of SED type. Note that the SWIRE library has only a 
few SEDs of starbursts. We think that a large variety of starburst 
SEDs, as in the SBURT model, is necessary to reproduce a good 
fraction of observed UV-mid-IR SEDs of IR-selected galaxies. 

In order to show what kind of sources are rejected in the 
SBURT fitting, we plot both accepted and rejected samples 
in a colour-colour plot of AG - S 7 vs N2 - N3 in Figure [4] 
Takagi et al. (2007a) use this colour-colour plot to discriminate 
AGNs from normal star-forming galaxies. AGNs identified with 
optical spectroscopy have red colours in both N2 - N3 and 
AG — SI. The accepted fits avoid this AGN colour region, be- 
cause the SBURT model does not have any AGN component. 
The accepted sample also avoids the clump of objects with blue 
colours at N2 - AG - 0.6 and AG - S7 -1.5, which are stars. 
The rest of the rejected sample have similar N2 - AG and AG -57 
colours to those of the accepted sample. These include quiescent 
spiral galaxies, which have systematically redder optical-near- 
IR colours than any of the SBURT models. As we discuss in 
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Section 5.1, limitations on the dust model may also cause rejec- 
tion of the best-fitting model. 

In order to derive any firm conclusions, photometric red- 
shifts of the good-fit sample should be accurate enough. Figure 
[3] shows the comparison of photometric redshifts to spectro- 
scopic redshifts. Solid circles and crosses represent galaxies 
with acceptable and rejected fits, respectively. Defining a catas- 
trophic error by > 0.2, we find no catastrophic errors and 

cr a.- = 0.034 for the good-fit sample, and 21/68 catastrophic 

a+z) 

errors for the bad-fit sample. Thus, the goodness of the SBURT 
fit is actually a good measure of the accuracy of photometric 
redshifts. 

In Figure [5] we also plot photometric redshifts de- 
rived from ground-based optical-NIR (BVRi'z' JK) bands 
only, using a widely-u sed photometric redshift code, Hyperz 
( Bolz onella et alj 20001). In Hyperz , we ad opt a population syn- 
thesis model of Bruzual & Chariot] d2003l) . which is distributed 
along with the cod e, and the extinction curve of starbursts 
dCalzetti et al.l 120001) for dust reddening. In order to compare 
the photometric redshifts with AKARI bands to those without 
AKARI, we show the good-fit sample (with SBURT fitting) only 
in Figure[5] We find one catastrophic error in ground-based pho- 
tometric redshifts, while this galaxy has a reasonable photomet- 
ric redshift with AKARI bands. Furthermore, ground-based pho- 
tometric redshifts seem to have systematic errors; i.e. at z ~ 0.5, 
most of photometric redshifts are lower than the spectroscopic 
redshifts. Such errors can be removed by using AKARI bands. 

In order to check the statistical reliability of our photo- 
metric redshifts, we us e a correlation betw een the NIR colour 
N2 - N3 and redshift (Tak agi et al.ll2007al) . With this test, we 
could detect systematic effects, if any, due to the usage of mid- 
IR fluxes for deriving photometric redshifts. In Figure [6] we 
show the correlations for both spectroscopic and photometric 
samples. For the photometric sample, we obtain the correlation 
of (N2 - N3) = (2.37 + 0.12)log(l + z) - 0.56, which is con- 
sistent with that of the spectroscopic sample, i.e. (N2 - N3) = 
(2.43 + 0.05) log(l +z) — 0.59. Therefore, we find no systematic 
errors in our photometric redshifts. 



4.2.2. PAH-selected galaxies 

We here show more details of SED fitting results specifically 
for PAH-selected galaxies. In Figure IA~T1 - IA.4I we show the 
best-fitting (both accepted and rejected) SED models for PAH- 
selected galaxies at z ~ 0.5 and 1 . We obtained acceptable SED 
fits for 22/38 and 10/18 of the PAH-selected galaxies at z ~ 0.5 
and 1, respectively. Despite the fact that PAH-selected galaxies 
are starbursts, the SBURT model does not provide acceptable 
fits for 40 % of PAH-selected galaxies, although the success rate 
of the fitting is better than that of whole sample. This topic is 
discussed in detail in the next section. 

We show the redshift distributions of PAH-selected galaxies 
in Figure [7] They are consistent with the expectation from the 
model templates, shown in Figure|2]and[3] There are no system- 
atic differences in photometric redshifts between the good- and 
bad-fit samples. This may indicate that the resulting large x 1 is 
not due to an incorrect redshift. The total-IR luminosity from the 
SED model is indicated in each panel of Figure lATTI — IA. 41 We 
note that PAH-selected galaxies at z ~ 1 include ULIRGs. This 
type of ULIRGs can be found only at high redshifts (see Section 

E2J. 



E 1 




Fig. 7. Redshift distribution of PAH-selected galaxies. Vertically 
and diagonally shaded histograms represent ll-to-7jum and 
15-to-9//m selected galaxies with good SED fits, respectively. 
Histograms with thin dashed and solid lines indicate the redshift 
distribution including the bad-fit sample. Spectroscopic redshifts 
are adopted if available. 
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Fig. 8. 8 pm luminosity and stellar mass as a function of redshift. 
All-band-detected sources and PAH-selected galaxies are repre- 
sented by small dots and open circles, respectively. We show 
galaxies with good SED fits only. 



5. Discussion 

5.1. Nature of PAH-selected galaxies 

In Figure [8] we show 8/im luminosities and stellar masses of 
PAH-selected galaxies, along with those of all-band-detected 
sources. The 8 /mi luminosity is derived from the observed flux 
in the IRC bands closest to the reset-frame 8 /mi and the lu- 
minosity distance from photometric redshifts. Some of PAH- 
selected galaxies at z ~ 1 have L% im > 10 n L G , which corre- 
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Fig. 9. Ratio of observed to model fluxes as a function of rest- 
frame wavelengths for PAH-selected galaxies at 0.4 < z < 
0.5. Open and solid circles indicate individual galaxies and the 
weighted average for each photometric band, assuming the mean 
redshift of 0.46 and 0.45 for good and bad fit sample, respec- 
tively. 



spond to the luminosity of ULIRGs for a typical starburst SED 
dWeedman & Houckll2008l) . 

Stellar masses are derived from the best-fitting SED model 
when it is accepted. The adopted initial mass function (IMF) is a 
top-heavy one with the power law index of x = 1.10 (the Salpeter 
IMF has x = 1.35), and the lower and u pper mass limits are 0.1 
and 60 M (see also lTakagi et al.ll2004 . For the Salpeter IMF, 
stellar mass increases by a factor of 2. Stellar masses thus de- 
rived should be regarded as a lower limit, since the contribution 
from evolved stellar populations is not included. Most of PAH- 
selected galaxies at z ~ 1 have stellar masses of > 3 x 1O 1O M . 
Even some fraction of PAH-selected galaxies at z ~ 0.5 have 
similar stellar masses. What are descendants of such massive 
PAH-selected, i.e. starburst galaxies? At redshift z ~ 2, mas- 
sive starburst galaxies ha ve been found ma inly as submillime- 
tre galaxies (SMGs; e.g. iBorvs et ail 120051) and they are good 
candidates of massive sph eroids (Blai n et al.l|2004t ISmail et al.l 
2004; iBouche et al1l2007t iTacconi et al.ll2008l) . However. PAH- 
selected galaxies are not as massive as SMGs. Their luminos- 
ity is typically an order of magnitude less than that of SMGs, 
corresponding to luminous infrared galaxies (LIRGfl). Recent 
high resolution images of LIRGs at z ~ 1 show that the 
majority of LIRGs ha ve disk morphology (Elbazetal. 2007; 
Melbourn e et al.ll2008l) . although the sample size is still small. 
Clearly, we need more comprehensive study to identify what 



triggers starburst activity of PAH-selected galaxies, in order to 
constrain their evolutionary path. 

What is the main cause of a large PAH-to-continuum flux 
ratio? As shown in Figure [2] and [3j the SBURT model predicts 
that galaxies with large PAH-to-continuum flux ratios are young 
starbursts with the starburst age of ^ 0.4 Gyr. For the whole sam- 
ple of PAH-selected galaxies with good SED fits, we derive the 
mean starburst age of 0.4 Gyr. We note that, in the SED fitting, 
the starburst age depends not only on th e PAH-to-continuum flux 
ratio, but also on the optical colours (Takagi et al. 1999). The 
lack of very young starbursts in the sample may be accounted 
for by a selection effect, i.e. such young systems are hard to be 
observed due to a short life time. 

Indeed, star-forming regions, i.e. young stellar systems 
in a galaxy have a large PAH-to-continuum flux ratio, as 
shown in Spitzer/IRAC ima ges of nearby galaxies in the 
literature. I Wang et ail (120041) study the Antennae galaxies 
(NGC4038/4039) and show the map of the 8.(fl-to-4.5yum flux 
ratio, which is close to the PAH-to-continuum flux ratio in our 
definition. This flux ratio is quite high only at the overlap re- 
gion of the two disks where the most of IR emission come from. 
They estimate that such active star-forming region comprise only 
~10% of the total stellar mass of the entire system; i.e. the ac- 
tivity is localized rather than global. On the other hand, PAH- 
selected galaxies should have global star-forming activity, since 
the flux ratio for the entire system is high. 

For fair discussion, we should also pay attention to the case 
with no acceptable fits and clarify the cause of poor fits. In 
Figure|9] we show the observed-to-model flux ratio as a function 
of rest-frame wavelength for PAH-selected galaxies at 0.4 < z < 

0. 5. All of the rejected models, but one in Figure[9]have a star- 
burst age of 0.6 Gyr, i.e. the oldest and reddest model in the tem- 
plate. This means that the model-fitting fails for galaxies with 
relatively red optical colours, but with strong PAH emission (see 
also Ohyama et al. in preparation). Since the SBURT model is a 
rather simple model of starbursts, which have a single starburst 
region and starburst stellar population only, it is not surprising to 
find its limitation. 

Figure [9] shows that a large discrepancy between observa- 
tions and the model can be found at the 5-7 /mi range with 
sharp rising trend towards longer wavelength. We can see the 
same trend for the good-fit sample as well, although it is less 
significant. This indicates that the dust model adopted in the 
SBURT may need to be modified. One possibility is the ioniza- 
tion of PAHs. The optical constant of PAHs in the SBURT model 
is for neutra l PAH s, which may not be true for every PAHs. 
iLi & Draind yOOl) suggest that ionized PAHs have systemati- 
cally higher absorption cross section at 4 - 10 /mi, compared to 
neutral PAHs. If some fraction of PAHs are ionized, the discrep- 
ancy found in the SED fitting may be recovered. 

Another possibility is that stars dominating optical light is 
not the same as those responsible for strong PAH emission. As 
starbursts age, dust heating would be mostly dominated by newly 
formed stars in molecular clouds, while most of optical light 
come from stars having already exited from molecular clouds; 

1. e. older stars are less attenuated and dominate optical light, 
while younger stars embedded in dense molecular clouds dom- 
inate infrared light. Such age-dependent dust attenuation would 
be more effective in the later stage of starbursts. This may ex- 
plain why we find larger discrepancy in optically red, probably 
old, starbursts. 



3 Galaxies with a total infrared (8 - 1000 /mi) luminosity of 
10"- 12 L o . 



4 Non-stellar emission, i.e. ~ 5 % of stellar emission is subtracted 
from the total flux 
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Above two possibilities could be discriminated with the far- 
IR photometry. The former, the case of ionized PAHs, would 
change only mid-IR luminosity and does not affect the far-IR lu- 
minosity. If the latter (age-dependent attenuation) is the case, the 
current model would underestimate not only mid-IR luminos- 
ity but also far-IR luminosity. Thus, expected PAH-to-total-IR 
luminosity depends on the scenario. Sensitive far-IR photome- 
try of red PAH-selected galaxies, e.g. with the Herschel Space 
Observatory, could have diagnostic power to identify the main 
cause of the discrepancy. 

5.2. PAH-to-total IR luminosity relation 

The PAH luminosity of local ULIRGs is less than that expected 
from the correlation of PAH and total IR luminosity for less lu- 
minous starburst galaxies. In Figure[10l we show this luminosity 
relation for galaxies at z < 0.2 with crosses. The PAH luminos- 
ity is represented with th e PAH 7.7/mi peak luminosity taken 
from Weedm an & Houckl (12008). The total infrared lumin osity 
is taken from the references in Weedman & Houckl J2008). and 
converted to the cosmology adopted here. The 7.7/mi peak lu- 
minosity of ULIRGs is systematically below the correlation for 
local starbursts. A most plausible explanation of this system- 
atic effect is an AGN contribution to the total IR luminosity. 
Ilmanishi etail (120071) claim that 30-50 % of optical non-Seyfert 
ULIRGs at z < 0.15 have dust-enshrouded AGN. Combined 
with optical Seyfert ULIRGs, their results i ndicate that ~ 50 % 
of ULIRGs at z < 0.15 harbor AGN (see also lGenzel et all 19981: 
iLutz et al.ll998l) . However, Figure[10]shows that none of z < 0.2 
ULIRGs are consistent with the pure starburst case. This may 
indicate another explanation for the lower 7.7 //m luminosity of 
ULIRGs. If the starburst region is heavily obscured, the absorp- 
tion of PAH emissions by silicate dust may not be negligible. 

High redshift ULIRGs may behave differently, owing to pos- 
sible evolutionary effects. We investigate the PAH-to-total IR 
luminosity relation of the all-band-detected sources, including 
PAH-selected galaxies. The 7.7/mi peak luminosities are esti- 
mated from the IRC photometry by the following meth od. For 
a local starburst spectral template of iBrandl et alj d2006l) . filter- 
convolved flux at the rest-frame 8 //m is approximately half of 
the 7.7 jum peak flux; 0.50 in the L15 filter at z = 1 and 0.42 for 
the 511 filter at z — 0.5 for example. Here we adopt 0.50 for 
the conversion factor, since z ~ 1 sources are our prime targets. 
This estimate would be accurate within a 30 % uncertainty, as far 
as the 8 /mi emission is dominated by starbursts. We note that 
intrinsic s pectral variation at 5 - 8 /mi is found to be small in 
starbursts dNardinietal] F2008). The rest-frame 8 //m luminosi- 
ties are estimated from the observed fluxes at the band clos- 
est to the rest-frame 8 /mi and the luminosity distance derived 
from the redshift estimates. We have not done any interpolation 
or ^-correction in this step, considering the uncertainty of the 
photometric redshift. The total IR luminosities are derived from 
the best-fitting SED model. Therefore, we consider the good-fit 
sample only. A well-known tight correlation between the global 
far-IR and radio emission allows us to check the reliability of 
estimated total-IR luminosity from the SED fitting. Among all- 
band-detected sources, we detected five galaxies in our 1 .4 GHz 
observation with WSRT (White et al. 2010 in prep), which cov- 
ers a part of the NEP-Deep field. Comparing the observed radio 
fluxes to those predicted from the best-fitting SBURT models 
based on the far-IR/radio correlation, we estimate that the de- 
rived total-IR luminosities have uncertainties of 50 %. Since the 
derivations of both the 7.7 fan and total IR luminosity assume 
that starbursts dominate the infrared emission, the results for 
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Fig. 11. Same as Figure|9]but for PAH-selected galaxies at 1.1 < 
z < 1 .4. Open and solid circles indicate individual galaxies and 
the average for each photometric band, assuming a mean redshift 
of 1.22. Dashed lines bracket the 3.3 /mi region. 



PAH-selected galaxies, the most plausible starburst candidates 
in our sample, should be most reliable. 

As shown in Figure[10] some PAH-selected galaxies at z ~ 1 
are ULIRGs. Their 7.7 /mi peak luminosities are larger than local 
ULIRGs with a similar luminosity, and consistent with the PAH- 
to-total IR luminosity ratio of less luminous starbursts. They 
have no local counterparts and we specifically call these PAH- 
selected ULIRGs "PAH-luminous galaxies". Total IR luminos- 
ity of individual galaxies can be found in the legend of Figure 
IA. 1 1 — |A~4~1 In our sample, we find no PAH-luminous galaxies at 
z ~ 0.5. 

Normal starbursts, which have unabsorbed PAH emission 
features and no AG N contribution, can reach h igher luminosity 
at higher redshifts. IWeedman & Houckl (1200 8) reach the same 
conc lusion by compiling various Spitzer IRS spectroscopic sam- 
ples. Rig bv et al.l ([2008) investigate the rest-frame 8 /mi lumi- 
nosity of galaxies at z ~ 2 as a function of total IR lum i nosity , 
and find a similar trend. Using Spitzer/IRS. lHuang et alj J2009) 
also report ULIRGs with high PAH-to-total IR luminosity ratio 
at z ~ 2. Our results show that there are extreme starbursts with 
ULIRG luminosity even at z ~ 1 . 

5.3. Photometric detection of PAH 3.3/im feature? 

For galaxies at z ~ 1, the PAH 3.3 //m feature falls into 57 
band. PAH-selected galaxies at z ~ 1 are good candidates 
for detecting t his feature, because of both redshift and PAH 
prominence. MagneUietaLl d2008l) made a similar attempt us- 
ing Spitzer IRAC photometry for galaxies at 0.6 < z < 0.9, and 
detected an excess over the stellar continuum, which may origi- 
nate from the PAH 3.3 fim feature. 

In Figure ITOl we can see that some objects, i.e. ID2711, 
2836, 2932, 3033, 4956, have blue SI - S9W colours, which 
might be caused by PAH 3.3 /mi feature. Surprisingly, this excess 
is notable even if compared with the model flux which already 
includes the contribution from 3.3 /mi feature. We show the com- 
parison of observed fluxes to model fluxes in Figure QT| Here 
we use the sample of PAH-selected galaxies at 1 . 1 < z < 1 .4 
with good SED fits. The average observed-to-model flux ratios 
range from 0.75 to 1.23, depending on the wavelength. At PAH- 
dominated wavelengths, i.e. > 7 /zm, the scatter is large and the 
ratios are marginally consistent with the model. At the rest-frame 
3 /mi band, including PAH 3.3//m feature, the flux ratio is sys- 
tematically high by a factor of 1 .23 on average. If this excess is 
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Fig. 10. 7.7/itn peak luminosity versus infrared luminosity. Solid and das hed lines represen t the correlation between these lumi- 
nosities found for less luminous starbur st galaxies and its 1 cr unce rtainty dHouck et alJl2007h . Legend describes symbols for local 
galaxies from Spitzer IRS sample from Weedm an & Houckl ((2008), all-band-detected sources, and PAH-selected galaxies at both 
Z ~ 0.5 and 1 . Solid circles indicate galaxies with z > 1 . Typical uncertainty of both luminosities is indicated at the bottom right 
corner. 



accounted for by a systematically strong 3.3 yum feature only, this 
feature must be unusually strong by a factor of 3 - 7, compared 
with the model. In the literature, no galaxies with such strong 
3.3 //m feature can be found to our knowledge. This anomaly 
could originate from the continuum emission, rather than PAH 
3.3//m feature. However, we stress that, for galaxies at z ~ 0.5, 
the SBURT model reproduces the rest-frame 1-3 fim contin- 
uum emission well, as shown in Figure [9] 

It may be possible that the model lacks additional very hot 
dust component which contribute to the rest frame 3-4 fim emis- 
sion. Such a near-IR excess with a colour temper ature of ~10 3 K 
is reported in stu dy of nearby galaxies (e.g. iLu et alj [2003; 
Flag ev et alj 1200 6). However, we note that galaxies with large 
3 fim flux tend to have small 4 /mi flux as well, compared with 
the model. If the additional component has a gray body emission, 
it is difficult to explain such a trend. Furthermore, PAH-selected 
galaxies at z ~ 0.5 show no systematic near-IR excess, compared 
to the model, as shown in Figure [9] 

We need more careful investigation of the origin of the flux 
anomaly at this wavelength range. In fact, the model does not 
take into account the absorption features due to Hydrogenated 
Amorphous Carbon (HAC) and/or H2Q. Unfortun ately, we have 
to wait for the l aunch of JWST (IClampinl I2008T) and SPICA 
dNaka gawa 2004) to obtain mid-IR spectra of these galaxies and 
reveal the origin of this flux anomaly. 



6. Summary 

Using a multi-wavelength mid-IR survey, the AKARI NEP- 
Deep, we have constructed a catalogue of all-band-detected 
sources. From this catalogue, we photometrically identified 
galaxies whose mid-IR emission is dominated by PAH emission. 
These PAH-selected galaxies have large PAH-to-continuum flux 
ratio, i.e. ll-to-7/im and 15-to-9,um flux ratios at z ~ 0.5 and 
1, respectively. PAH-selected galaxies are the best candidates 
of starburst galaxies at these redshifts. Some of PAH-selected 
galaxies have stellar masses of > 3 x 1O 1() M , and could be 
progenitors of present-day massive galaxies, which are mostly 
spheroids. 

An evolutionary SED model of starbursts, SBURT, repro- 
duces observed optical-to-mid-IR SEDs of more than half of 
PAH-selected galaxies. According to the SED model, the PAH- 
to-continuum flux ratio of > 8 can be explained with the starburst 
age of 0.4 Gyr or younger. The average starburst age from the 
SED fitting is 0.4 Gyr, which is marginally consistent with the 
PAH-to-continuum flux ratio. The lack of young starbursts may 
be the selection effect due to a short life time of young starbursts. 
On the other hand, SBURT has a difficulty to reproduce large 
PAH fluxes of optically red PAH-selected galaxies, which are 
probably evolved starbursts. This may require the SED model to 
include the age-dependent extinction, and/or to adopt improved 
optical properties of PAHs. 
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At z ~ 1, the infrared luminosity of some PAH-selected 
galaxies corresponds to that of ULIRGs, and we call them PAH- 
luminous galaxies. They have large PAH luminosity, compared 
with local ULIRGs, but the PAH-to-total IR luminosity ratio is 
comparable to that of less luminous starbursts. PAH-luminous 
galaxies seem to be a unique galaxy population at high redshifts. 
The number density of PAH-luminous galaxies would be given 
elsewhere, using whole area of the NEP-Deep survey. 

There is a hint that the PAH 3.3 //m feature is quite strong 
in PAH-luminous galaxies by a factor of 3 - 7, but this needs to 
be confirmed with next generation infrared telescopes, such as 
JWST and SPICA, with infrared spectroscopy. 
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Appendix A: SED fitting results 

Here we show the results of the SED fitting analysis of PAH- 
selected galaxies in Figure |A~T1 - IA.4I and present the AKARI 
photometry in Table A. 1 . One source (ID4475) has problematic 
photometry in the S-cam images, due to a nearby bright source, 
and therefore is not included in the figure. 
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A.2. Same as Figure IATT1 but for the bad-fit sample. 
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Fig. A.3. Same as Figure IATT1 but for PAH-selected galaxies at z ~ 1 with acceptable fits. 
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A.4. Same as Figure IA31 but for the bad-fit sample. 



Table A.l. AKARI coordinates and flux densities of PAH-selexted galaxies in the S-cam field 
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Notes. — Galaxies with ambiguous optical identification are excluded. Column (1): AKARI mid-IR source ID. Columns (2): Source Name. Columns (3) and (4): AKARI W2-band J2000.0 RA and Dec. Columns (5) through 
(22): IRC flux densities and errors. 



Table A.l. - continued 



ID 


Name 


RA 


Dec 


2.4/im 


3.2 jum 


4.1 ftm 


l.Qfim 


9.0/im 


11 fim 


15//m 


18^m 


24 fim 








[J2000] 


Jv 


A f 


Jv 


A f 


Jv 


A f 


Jv 


A f 


Jv 


A f 


Jv 


A f 


Jv 


A f 
&Jv 


Jv 


A f 


Jv 


A f 












fiJy 




/Jy 


/Jy 




/Jy 


/Jy 


/Jy 


/Jy 


/Jy 


/Jy 


/Jy 


/Jy 


/Jy 


/jjy 


/ijy 


fily 






(1) 


(2) 


(3) 


(4) 


(->) 


\P) 


f"7~l 


\p) 




n c\\ 

[IVi) 


(\ 1 ^ 


("1 7\ 


^13) 






(\ fC\ 
Uoj 


n 7"* 


Uo) 


/I Q\ 

uy; 


(ZUj 


{LI) 






3110 


J175331.90+663231.9 


268.38295 


66.54220 


OJ. 1 


Z.3 


oi. y 


1 o 

i.y 


/in 7 
4U. / 


1 ^ 
l.J 


11/1 


o.U 


A 7 7 

03.3 


1 n 7 
1U. / 


77n n 
Z3U.U 


i n n 

iy.u 


i /i n fi 

i4y.u 


7n /i 
ZU.4 


1 Q/i n 
104.U 


77 Q 

ZZ.o 


7 i c n 
Zlo.U 


40. / 




3387 


J175345. 13+663405.1 


268.43805 


66.56809 


JZ.D 


7 ^ 


AA A 
44.4 


7 7 
L.L 


71 1 

31 .3 


1 7 


Q 7 

y. / 


y ^ 
o.j 


y3.o 


11 Q 

1 1.0 


7<n n 


in /i 

iy.4 


1 7/1 n 
1 /4.U 


7/1 n 
Z4.U 


1 7/i n 
134.U 


7n 7 
ZU.3 


77 7 
33. / 


4 J. 4 




3953 


J175456.90+664818.7 


268.73711 


66.80520 


G'J 7 

yj.Z 


L. 1 


Q/1 Q 


7 7 
2.3 


A7 1 
02. 1 


1 C 

1.8 


1 1 c 

1 l.J 


Q 7 

o.Z 


1 on n 


1 n a 
1U.4 


71 1 n 
31 l.U 


1/11 

14.1 


7nc n 
jUj.U 


77 7 
Z /.3 


7Q/i n 

zy4.u 


71 O 

zi.y 


7^A n 


DU.3 




4359 


J175629.33+664530.6 


269.12221 


66.75850 


<7 A 


7 n 
L. 1 


04. U 


7 7 
L.L 


AO 7 

4y. / 


1 7 
1. / 


1 /.j 


5. / 


7/1 7 

Z4.3 


Q C 

O.J 


men 


1 7 < 
12. j 


7Q7 n 

joj.U 


33. J 


iqc n 
Ioj.U 


7n < 
3U.J 


/1Q7 n 
45Z.U 


1 1 i n 
1 1 l.U 




4438 


J175319.00+664419.7 


268.32918 


66.73881 


1 lo.O 


7 n 

2.y 


i c\i n 
103.0 


7 1 


(TO A 


1 "7 
1. / 


on 7 
ZU.Z 


7 O 
1.1 


i cn n 
IjU.U 


1 n i 
10.3 


OKA 

zl j.O 


1 a a 
13.4 


Zlfi.O 


on n 

zo.y 


7yi c n 
24 j .0 


m i 

iy.i 


7 no n 
30fi.O 


/in i 
49.1 




4454 


J175335.34+664522.4 


268.39727 


66.75624 




2. / 


A7 1 
/. 1 


7 1 
2. 1 


zl7 Q 

4 /.y 


1 7 
1. / 


70. 8 
2U.0 


7 7 


SO 8 

sy.o 


1 1 1 
ii.i 


70.7 n 


1/1 O 
14. U 


7^1 ( \ 
Zj 1 .U 


71 Q 

zi.y 


77H n 
2 /U.U 


1 Q 7 

ly.z 


7AS n 
ZOo.U 


AS 1 
40. 1 


H 


4653 


J175645.64+664334.9 


269.19020 


66.72637 


82.1 


3.2 


74.2 


2.4 


51.6 


2.2 


37.6 


8.1 


296.0 


10.1 


463.0 


12.9 


333.0 


27.8 


347.0 


23.3 


465.0 


64.4 




4802 


J175328.00+664228.0 


268.36667 


66.70779 


66.5 


2.3 


67.2 


2.2 


43.9 


2.0 


8.9 


6.6 


60.5 


9.2 


200.0 


11.2 


245.0 


19.8 


132.0 


17.7 


143.0 


44.7 




5339 


J175651. 89+664003.0 


269.21623 


66.66751 


50.0 


2.4 


45.8 


2.2 


28.8 


2.1 


7.8 


9.7 


60.2 


10.5 


217.0 


13.2 


209.0 


32.0 


216.0 


29.5 


247.0 


62.9 


"0. 


5424 


J175534.22+663918.4 


268.89262 


66.65512 


53.7 


2.0 


65.1 


1.9 


50.2 


1.4 


42.7 


11.4 


136.0 


12.7 


538.0 


19.1 


487.0 


35.9 


522.0 


47.6 


1010.0 


96.6 


S- 


5776 


J175328.61+663701.4 


268.36922 


66.61707 


53.6 


2.0 


48.8 


1.8 


31.0 


1.4 


3.2 


7.0 


75.2 


10.8 


162.0 


16.3 


200.0 


19.8 


195.0 


18.2 


250.0 


47.6 




6363 


J175701.56+665331.4 


269.25654 


66.89206 


54.1 


2.4 


61.1 


2.3 


39.0 


1.7 


4.4 


7.4 


86.4 


9.8 


197.0 


13.2 


270.0 


25.6 


278.0 


30.6 


507.0 


86.5 




6393 


J175402.88+665222.2 


268.51203 


66.87284 


257.0 


4.5 


229.0 


3.3 


174.0 


2.8 


176.0 


10.2 


828.0 


11.9 


1420.0 


17.9 


1190.0 


29.8 


1400.0 


28.7 


1780.0 


62.0 


O, 


6427 


J175430.88+665249.1 


268.62869 


66.88033 


141.0 


3.1 


132.0 


2.7 


89.1 


2.1 


65.5 


7.9 


342.0 


9.4 


630.0 


14.9 


636.0 


26.4 


537.0 


22.1 


653.0 


50.6 


■5 


6449 


J175542.96+665248.9 


268.92901 


66.88027 


37.2 


1.9 


39.3 


1.7 


27.6 


1.4 


7.5 


7.8 


72.5 


8.4 


213.0 


12.5 


232.0 


28.0 


274.0 


23.5 


426.0 


57.9 





Notes.— Column (1): AKARI mid-IR source ID. Columns (2): Source Name. Columns (3) and (4): AKARI /V2-band J2000.0 RA and Dec. Columns (5) through (22): IRC flux densities and errors. 
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